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Abstract 

Momentum distributions of particles from nuclear break-up of fast three- 
body halos are calculated consistently, and applied to 11 Li. The same two- 
body interactions between the three particles are used to calculate the ground 
state structure and the final state of the reaction processes. We reproduce the 
available momentum distributions from n Li fragmentation, together with the 
size and energy of 11 Li, with a neutron-core relative state containing a p-state 
admixture of 20%-30%. The available fragmentation data strongly suggest an 
s-state in 10 Li at about 50 keV, and indicate a p-state around 500 keV. 
PACS numbers: 25.60.+V, 21.45. +v, 21.60. Gx, 27.20. +n 
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Introduction. The existence of halo nuclei is by now well established as weakly bound 
and spatially extended systems [[J . The dominant features of these unusual nuclei are most 
easily described by few-body models. The many degrees of freedom are then divided into the 
approximately frozen (core) and the active (halo) degrees of freedom. Substantial efforts have 
been allocated to investigations of three-body halos (one core and two halo particles) and 
in particular to bound three-body systems where all two-particle subsystems are unbound 
||. The nuclear prototypes of these so-called Borromean systems are 6 He ( 4 He+n+n) and 
n Li ( 9 Li+n+n), both thoroughly discussed in a general framework in ||. 

The most detailed source of experimental information about the structure of these nuclei 
is the momentum distributions of the "particles" resulting from fragmentation reactions 
0-0]. One major problem in the interpretation of such measurements is the inherent mixture 
of effects from the original structure and the reaction mechanism || . Both should therefore 
be consistently incorporated in model calculations. 

The simplest model assumes the sudden approximation, where one of the particles in- 
stantaneously is removed from the three-body system while the other two particles remain 
completely undisturbed. Clearly this can only be justified for reaction times much shorter 
than the characteristic time for the motion of the three particles in the system. Since halo 
nuclei are weakly bound and the beam energy is very high this requirement is exceedingly 
well fulfilled. The observed momentum distributions therefore seems to provide direct in- 
formation about the three-body wave function |J. However, the interactions between the 
remaining two particles are often essential especially when low lying resonances are present 
p jl0jll|] . In principle, processes in which the three particles move along together after the 



fragmentation are also possible. However they occur through Coulomb dissociation, requir- 
ing a heavy target or by using of a low energy beam. Since we consider a light target and a 
high energy beam this kind of reactions is then out of our model, and they are not considered 
in this letter. 

The final state two-body interactions (FSI) are active in the last part of the reaction pro- 
cess and they determine simultaneously the three-body structure of the initial halo nucleus. 



A consistent treatment is necessary to allow reliable interpretations of this major source 
of detailed experimental information. In this letter we report on such model calculations, 
where we combine an accurate three-body description with an equally accurate computation 
of momentum distributions. After a general discussion we shall present detailed calculations 
of momentum distributions of particles from fragmentation reactions of n Li when FSI are 
included. 

Method. After the collision in which one particle is suddenly removed, the probability 
of finding the remaining two particles with relative momentum and total momentum h y 
relative to the center of mass of the three-body system is proportional to the overlap 

^(k x ,k y ) = (e ik ^V k - x |vI/(x,y)), (1) 

where \&(x, y) is the three-body wave function. The coordinates x and y are the usual 
Jacobi coordinates @J3|] where x is drawn between the two particles surviving after the 
fragmentation. Without FSI the momentum distribution or the differential cross section is 
then proportional to the square of the Fourier transform of the three-body wave function. 

Inclusion of FSI now amounts to substituting the plane wave e lkl ' x in eq. QT|) by the 
appropriate distorted two-body wave function w. The momentum distribution is then given 
by 
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where J is the total spin of the halo nucleus, and s x and a x are the spin of the two-body 
final state and its projection. The summations in eq.(|2|) arise from the average over initial 
states (M) and the sum over final states (s x and a x ). 

In our calculation the three-body wave function is obtained by solving the Faddeev 
equations in coordinate space, where the nucleon-nucleon potential is fitted to low energy s 
and p-wave nucleon-nucleon scattering data, and the neutron-core potential is adjusted to 
give the proper binding energy and mean square radius of the three-body system [|T2] , |I3]] . 

The partial wave expansion of the two-body final state wave function is written as ]14| 



w SCT (k,x) = 




where the radial functions u 3 is (k,x) are obtained by solving the Schrodinger equation with 
the appropriate two-body potential. Finally, we calculate the expression in eq.(|2|) and sub- 
sequently integrate over the unobserved quantities to obtain the measured momentum dis- 
tributions. The details of the formalism will be presented elsewhere. 

General results. We reveal the general results of the procedure described above by using 
the three-body halo n Li ( 9 Li+n+n) as an example, and neglecting the spin of the core. The 
neutron-neutron interaction is fitted as previously indicated, and the neutron-core s and 
j9-wave potentials are varied to study the dependence of the momentum distribution on the 
positions of the resonances in 10 Li, as well as the effect of the inclusion of FSI. The binding 
energy and mean square radius of 11 Li is kept fixed to the experimental values. 

In the upper part of Fig.l we show the resulting neutron momentum distributions after 
one of the neutrons is removed by the target for three different energies of the sy2 virtual 
state (100, 200 and 300 keV). The upper set of curves is obtained without FSI. They show 
a little variation, since they depend mainly upon the binding energy and r.m.s. radius of 
11 Li. The final state interactions substantially reduce the widths, the stronger the lower 
the virtual state is. In table I we summarize the full width at half maximum (FWHM) for 
the curves in Fig.l. Here we also give the FWHM of the core one dimensional momentum 
distributions, which are much less influenced by FSI due to the larger mass. 

The lower part of Fig.l shows the same kind of calculation for core break-up reactions, 
where the core is destroyed during the interaction, keeping the validity of the sudden ap- 
proximation. The interest in this kind of reactions is that now the final state interaction 
(neutron-neutron interaction) is well known, reducing the uncertainties involved in the cal- 
culation. We first observe that variation of the position of the S\ji virtual state only pro- 
duces a small change in the momentum distributions. That is expected, since the final state 
neutron-neutron interaction is independent of the properties of 10 Li (n+core). However when 
comparing with the results without FSI, as shown by the upper set of curves, we still observe 
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a noticeable variation (note that the curves without FSI are identical to the corresponding 
curves in the upper part of the figure). Now the FWHM of the distribution is around 40 
MeV/c, broader than the distributions obtained when one of the neutrons is removed, but 
narrower than the calculation without FSI. As seen from table I, the agreement with the 
experimental data is quite good, especially for low energies of the s virtual state. This is in 
complete agreement with the low lying s virtual state (~ 50 keV) recently suggested in f| 
along with a p-state resonance near 0.5 MeV. 

It should be mentioned that a narrow momentum distribution can also be obtained with 
a low lying p-resonance. However in this case the lowest s-state must be shifted towards 
much higher energies to reproduce the correct binding energy of 11 Li. 

Realistic calculations for 11 Li. Up to now, we have neglected the spin dependence of the 
neutron-core interaction or equivalently assumed that the spin of the core is zero. Of course, 
this is not realistic (the spin of 9 Li is 3/2), and the neutron-core interaction should include 
a spin dependence, splitting the two possible s-states of 10 Li with total angular momentum 
1 and 2. To do this, we have in the neutron-core potential included one term proportional 
to s„ • s c , where s„ is the spin of the neutron and s c the spin of the core [13]. For simplicity, 
the spin splitting term has been introduced only in the s wave. 

It is now possible to place an s-state at 50 keV and the lowest p-resonance at 500 keV, 
as suggested by the analyses in ||. Simultaneously we are able to vary the content of 
s and p-waves in the neutron-core subsystem of 11 Li. Such realistic calculations can be 
compared directly with the experimental data. We then place the lowest p-resonance at 
500 keV, and the lowest s-state with total angular momentum 2 at 50 keV. Then we use 
the spin-orbit potential in the neutron-core interaction to vary the total p-state content in 
the n Li wave function, and finally we adjust the energy of the second s-state (with total 
angular momentum 1) to recover the correct binding energy and mean square radius in the 
n Li. This procedure determines completely the low lying resonance structure of 10 Li. Very 
similar results are obtained by placing the s-state with angular momentum 1 at 50 keV, and 
instead adjusting the energy of the s-state with angular momentum 2. 



In Fig. 2 we show the one dimensional transverse momentum distribution of 9 Li fragments. 
The experimental data [|7| correspond to a reaction at 280 MeV/u in a C target. The n Li 
wave function contains around 26% of p-wave in the neutron-core subsystem, but the result 
is almost independent of the p-state admixture. The results of the computation with (solid 
line) and without (dashed line) FSI are shown in the figure. The curves have been convoluted 
with the experimental beam profile [0]. As expected, the effect of FSI is rather small in the 
central region. For momenta larger than around 75 MeV/c the neutrons are inside the core, 
and our three-body model is not valid anymore. Also, in the high momentum region other 
effects, like diffraction processes, should be considered. The good agreement in the central 
region shows that our n Li wave function is accurate enough to describe the process. 

To evaluate the effect of FSI on the momentum distributions we investigate neutron 
distributions, much more sensitive to them. The results are shown in Fig. 3 for the two 
dimensional neutron momentum distribution, both for neutron removal process (upper part) 
and for core break-up reactions (lower part). The appropriate FSI are included, i.e. neutron- 
core in the first case and neutron-neutron in the second case. In both cases the experimental 
data correspond to reactions of 11 Li at 280 MeV/u with a C target Several 11 Li wave 

functions with different content of s and p-waves in the neutron-core channel have been used. 
In particular, going from the narrower to the broader distributions, the curves correspond 
to calculations with 4, 18, and 35% of p-wave. In both parts of the figure, FSI are essential 
to recover the observed behaviour of the distributions. 

The core break-up process (lower part) is almost insensitive to the structure of n Li and 
the peculiarities of the neutron-core interaction. This is due to the fact that the final state 
interactions between the two neutrons are independent of the neutron-core potential. The 
nice agreement with the experiment strongly supports the method and the model. In the 
upper part of Fig. 3 the two dimensional neutron momentum distribution for neutron removal 
reactions are shown. Now we see that a n Li wave function with a small content of p-wave 
clearly underestimates the width of the momentum distribution, while a wave function with 
more than 35% of p-wave in the neutron- 9 Li subsystem overestimate the width. The best 
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results are obtained when a p-wave admixture of around 26% is used. The data in the upper 
part of Fig.3 also support the choice of a 500 keV p-state, since this energy corresponds to 
the bump observed around p r ~ 30 MeV/c. 

Summary and conclusions. We have calculated momentum distributions of fragments 
from high-energy nuclear break-up reactions of three-body Borromean halo nuclei. The 
two-body interaction in the three-body description of the halo nucleus is identical to the 
final state interaction between the two particles remaining after the target has removed the 
third particle. We maintain this consistency in the description. The final state interaction 
may significantly reduce the width of the momentum distribution of the light particle. The 
size of this reduction strongly depends on the resonance structure of the remaining two-body 
system. 

After extraction of the general properties, we applied the model to fragmentation of the 
three-body halo nucleus 11 Li. First, we compared the core-momentum distribution after 
neutron removal with experimental data. The results are only marginally sensitive to final 
state interactions and the sudden approximation only recovers the ground state wave func- 
tion in momentum space. The good agreement with measurements therefore demonstrates 
the validity of our 11 Li three-body model and the reaction mechanism assumed. 

Secondly, we computed two dimensional neutron distributions after core break-up. They 
are rather insensitive to the characteristics of the neutron-core interaction which only enters 
through the properties of the ground state wave function. However the final state neutron- 
neutron interaction is essential to reproduce the measured data. This demonstrates the 
validity of our method of including final state interactions. 

Thirdly, we computed two dimensional neutron distributions after one neutron removal. 
They depend crucially on the neutron-core interaction which now enters both through the 
properties of the ground state wave function and through the final state interaction. The 
structure of n Li must then contain a p-state admixture of about 26% in the relative neutron- 
core system. The same interaction determines the structure of 10 Li, where a p-state resonance 
at about 500 keV produces the small shoulder at about 30 MeV/c in the experimental 
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distribution. The statistical average of the two spin splitted s virtual states must then be at 
about 700 keV in order to reproduce the binding energy of 11 Li. Furthermore, the narrow 
momentum distribution can only be reproduced with the lowest s virtual state in 10 Li at 
about 50 keV. 

Acknowledgments One of us (E.G.) acknowledges support from the European Union 
through the Human Capital and Mobility program contract nr. ERBCHBGCT930320. 



8 



REFERENCES 

[1] P.G. Hansen, A.S. Jensen and B. Jonson, to be published in Ann. Rev. Nucl. Part. Sci. 
45 (1995). 

[2] D.V. Fedorov, A.S. Jensen, and K. Riisager, Phys. Rev. C49, 201 (1994). 

[3] M.V. Zhukov, B.V. Danilin, D.V. Fedorov, J.M. Bang, I.S. Thompson, and J.S. Vaagen, 
Phys. Rep. 231, 151 (1993). 

[4] R. Anne et al, Phys. Lett. B250, 19 (1990). 

[5] M. Zinser et al, Phys. Rev. Lett. 75, 1719 (1995). 

[6] T. Nilsson et al, Europhys. Lett. 30, 19 (1995). 

[7] F. Humbert et al, Phys. Lett. B347, 198 (1995). 

[8] C.A. Bertulani, L.F. Canto, and M.S. Hussein, Phys. Rep. 226, 281 (1993). 

[9] I.S. Thompson and M.V. Zhukov, Phys. Rev. C49, 1904 (1994). 

[10] F. Barranco, E. Vigezzi, and R.A. Broglia, Phys. Lett. B319 387 (1993) 

[11] A. A. Korsheninnikov and T. Kobayashi, Nucl. Phys. A567, 97 (1994). 

[12] D.V. Fedorov, A.S. Jensen, and K. Riisager, Phys. Rev. C50, 2372 (1994). 

[13] D.V. Fedorov, E. Garrido , and A.S. Jensen, Phys. Rev. C51, 3052 (1995). 

[14] R.G. Newton, Scattering Theory of Waves and Particles, (Springer- Verlag, N.Y., 1982), 
Second Edition, p. 444. 



9 



TABLES 







neutron 






core 


no FSI 


FSl( n " c ) 


FS j(n-n) 


no FSI 


FSl( n - c ) 


E Sl/2 = 100 keV 


51.1 


28.3 


38.4 


68.5 


55.1 


E Sl/2 = 200 keV 


53.2 


33.9 


39.0 


69.8 


58.0 


E Sl/2 = 300 keV 


55.8 


39.9 


39.7 


71.4 


60.4 


Experimental 




25-35^ 


43 ± 3( b ) 




49 ± 3( c ) 



TABLE I. Comparison of the full width at half maximum in MeV/c of the longitudinal momen- 
tum distributions from n Li fragmentation for different positions of the S1/2 virtual state. Columns 
2 to 4 refer to neutron momentum distributions, while 5 and 6 refer to core momentum distribu- 
tions, (n-c) FSI between 9 Li and the neutron have been included (neutron removal process), (n-n) 
FSI between the two neutrons have been included (core break-up process), (a) See, for instance, 
|]. (b) Data from @. (c) Data from §. 
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Figure Captions 



Fig.l.- Upper part: The neutron momentum distributions from neutron removal process 
in n Li fragmentation. p x is the x-cartesian component of the momentum of the neutron 
relative to the center of mass of the three-body system. The syz virtual state is placed at 
100 keV (solid line), 200 keV (dashed line), and 300 keV (long dashed line), respectively. 
The upper set of curves corresponds to calculations without FSI, while FSI are included in 
the lower set. Lower part: The same as before for core break-up reactions. 

Fig. 2.- Transverse core momentum distribution after neutron removal from n Li. p x is 
the x-cartesian component of the momentum of the core relative to the center of mass of 
the three-body system. The solid (dashed) line is the calculation with (without) FSI. Data 
from 0. 

Fig. 3.- Radial neutron momentum distribution for neutron removal (upper part) and 
core break-up (lower part) reactions of n Li. p r = <Jp x + p^ is the cylindrical radial com- 
ponent of the momentum of the neutron relative to the center of mass of the three-body 
system. Results with (solid line) and without (dashed line) FSI are presented. The p wave 
content in the neutron- 9 Li subsystem (from narrower to broader distributions) is 4, 18, and 
35%. Experimental data are taken from || and ||. 
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